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1. INTRODUCTION 

The key to the success of a business is to be competitive. Now, to meet customer images over time, 
ensuring high quality products and services is not a difference, but a basic form. In this case, in general, 
with the best or near-optimal processes reducing costs, improving slowing and economic efficiency, 
particularly due to better use of product quality, resources, and better service levels, can be the shortest time 
and inventory levels available very complicated theme. As production plans increase, the complexity and 
hassle will need to create a master plan, grow faster, especially in the case of limited resources, as is the case 
in most companies [1]. In fact, Production scheduling problems are proved to be NP-hard sorts of problems 
and are not easily solved for large sizes [2]. Fortunately, new ideas and methods of improving production 
management, such as the use of heuristics AI, are frequently being proposed by researchers. In order to solve 
NP-hard issues, use of meta-heuristic methods requires both effective and efficient optimization techniques. 
A review of existing literature reveals that much work has yet to be reported in the application of meta- 
heuristic techniques to solve MPS problems [3]. Research in this field is continually being conducted, 
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with widely used nature-inspired meta-heuristic optimization techniques proving superior to traditional 
techniques. Although evolutionary computation methods provide solutions combining computational 
efficiency and good performance, evolutionary computational research has been criticized for considering 
artificial test problem scenarios that are much less complex than real-world manufacturing cases [4]. 

Multi-objective optimization is the concurrent technique that defines more than one contrasting goal 
under specific restrictions [5]. Use (MOEA) to solve multifunctional optimization issues. Some MOEA are 
population-based optimization, where repetition is performed on a group of solutions (population) and in 
each generation, multiple solution will return. For the following reasons MOEA is more popular to solve 
multi-objective optimization issues: (i) ease of implementation, (ii) robust and flexible, (iii) less chances of 
falling into local minimum, (iv) it gives the number of optimal solutions [6]. 

The researchers used a variety of advanced optimization techniques; Vieira et al [7] applied 
simulated annealing, Soares et al [1] proposed a different framework of genetic algorithms, Vieira [8] has 
compared simulated annealing and genetic algorithms for master production scheduling problems and 
Radhika et al [3, 9] applied differential evolution. The goal consideration was to minimize inventory level, 
maximum service level and minimize each of overtime and inventory level below safety stock. 

Wu et al [11] suggested BP neural network to solve the demand predicate for master production 
scheduling. The goal is considered to maximize device usage and minimize loss for penalties for production 
inventory. 

Bakar, et al [12] developed solutions on Master Production Schedule (MPS) with the help of an 
algorithm that is based on the installation of two birds swarm algorithms namely Bat Algorithm (BA) and 
Gravitational Search Algorithm (GSA). A new multi-objective evolutionary algorithm Bat Algorithm and 
Gravitational Search Algorithm MOGSABAT was also proposed that considered the configuration 
performance of the system parameters. Therefore, solving the issue of the MPS using intelligent algorithms 
were carried out by a specific number of researchers. 

A memetic algorithm (MA) has been proposed to solve the problem of multi-objective master 
production scheduling. Six scenarios have taken from previous studies and the result obtained are compared 
with two meta-heuristic algorithm simulated annealing and genetic algorithm. 

The paper’s framework is structured as follows: Section 2 illustrate a brief revision of production 
planning problem, and the proposed algorithm has been explained in detail in Section 3. Section 4 shows the 
computational results, and finally, the outcome of this research has been provided in Section 5. 


2. OVERVIEW OF PRODUCTION PLANNING OPTIMIZATION AND MEMETIC 
ALGORITHM 
This section explains some basic concepts of production planning, particularly, master production 
scheduling and memetic algorithms. 


2.1. Master Production Scheduling 

As per the American Production and Inventory Control Society (APICS), a master production plan 
is a statement of what the company needs to produce, which in turn becomes a sequence of scheduling 
choices that guide the Material Requirements Planning (MRP) scheme. It represents what the company 
intends to produce and is expressed in configuration, quantities and specific dates. The master plan is not a 
sales forecast, which represents a demand declaration. It should take into consideration many others such as 
the requests, material availability, pending orders, projected ending inventory levels, capacity availability, 
managerial policies and goals, among others. 

MPS problems usually bring about conflicting objectives, such as minimizing of inventory, 
maximizing of service levels, and set-up times. For these reasons, it is suggested that heuristics and meta- 
heuristics be used to solve these kinds of issues. Several artificial intelligence meta-heuristics were used for 
optimization, including tabu search, genetic algorithms, simulated annealing, ant colony and beam search [7]. 


2.2. Mathematical Model of MPS 
The mathematical model of MPS problem can be represented as a mixed integer program as 
follows [1]: 


K P 
Lai Lp=1 El kp 
TH 


Min EI = (1) 


K P 
Le=1 Lp=1RNMkp 
TH 


Min RNM = (2) 
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K P 
Èk=1 Èp=1 BSSkp 


Min BSS = e 
TH 

Min OC = YR, &b=1 OCrp j 

Subject to: 

OH, if(p = 1) 
Bly =| ; 
kp Elkp-1) if (p > 1) i 
Elxp = max |o, ((MPST yp + Blyp) g GRxp)| m 
MPST yy = X3-1 MPSkpr k 
MPSppr = BNupy * BSicpr (8) 
RNMyp = max |0, (GRyp — (MPST yp + Blkp) )| (9) 
BSSip = max{[0, (SSkp — Elxp)| ve 
CUH,, = SK, OES) ae 
rp Kak URky 

CUH,p < ACrp i 
OC,» = max|0, (CUH,;» ~ AC,p)] oe 


2.3. Memetic Algorithm 

In this formulation, each character expresses the value as follows. K: Total quantities of separate 
products, R: Total quantity of separate productive resources, P: Total number of planning periods, TH: Total 
horizon of planning, E/kp: Ending inventory level generated for product k at period p, RNMkp: 
Requirements not met for product k at period p, BSSkp:Quantity bellow safety inventory level for product k 
at period p, OCrp: Over capacity needed at resource r at period p, Blkp: Initial inventory level of the 
product k at period p, OHk: Initial storage available (on-hand), in the first planning period, GRkp: Gross 
requirement for product k at period p, BSkp: Standard lot size for product k at period p, Nkp: Net 
requirement for product k at period p, considering infinity capacity, SSkp: Safety inventory level of product 
k at period p, URkr: Production rate for product k at resource r (units per hour), ACrp: Available capacity, 
in hours, at resource p at period p, BNkpr: Quantity of standard lot sizes needed for the production of the 
product k at resource r, at period p (number of lots), MPSkpr: Total quantity to be manufactured of the 
product k at resource r at period p, MPSTxp: Total quantity to be manufactured of the product k at period p 
(considering all available resources), CUHrp: Capacity used from the resource r at period p. 

Headings In the mid of 1960s and 1970s, various researchers laid the foundations of the 
evolutionary algorithms (EAs) we now know. Since this period, and despite some difficult beginnings, 
most researchers have been interest in searching or optimizing existential and theoretical perspectives to 
understand and accept the existence and authenticity of these technologies. The same is true for other related 
techniques, for instance simulated annealing (SA), tabu search (TS), and the like. The name metaheuristics is 
used to collectively term these specific technologies. At the end of the 1980s, the term "memetic algorithms" 
(MAs) was given to refer to meta-heuristic family, and they tried various concepts of close-up disintegration- 
in those days-families like EAs and SA. The word "memetic” comes from the term "meme" written by Dougs 
in 1976 for the analogy of genes in the context of cultural evolution [14]. 

The memetic algorithm was named by Moscato [15] after a series of stochastic global search 
techniques that were incorporated into the framework of evolutionary algorithms in a wide range of lines, 
"problem-specific local search heuristics" and multi-agent systems [16]. The web treatment method using 
local search provides several names in the research paper. The most meme algorithm is the mixture of local 
search heuristic and genetic algorithm. This method is also called genetic local search, hybrid genetic 
algorithm, hybrid evolution algorithm. And other names. This hybrid approach has been applied to a large 
number of optimization problems with only great success [17]. For example, the memetic algorithm for 
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multi-objective optimization has been used to address home healthcare [18] and Constrained Weapons 
assignment (WTA) problem [19]. Different to traditional evolutionary computation (EC) methods, 
MAS fundamentally concern the exploding of all existing knowledge about the problem being studied. 
The problem domain of problem solving is not an optional mechanism, but the basic characteristics of MAs's 
features [20]. 

MAs are like a EAs theory based on population. This implies that the algorithm holds a number of 
solutions for the problem at hand, in particular a pool of multiple solutions at the same time. These solutions 
are inspired by EA terminology, following the nature on which these technologies are based. In the context of 
MAs the two terms will be used interchangeably. Each individual provides a prudent solution to the problem 
being solved. These solutions are subject to the process of reciprocal cooperation and competition, 
which seems to make the behavioral patterns of similar individuals clearer at this point. The first thing to 
consider is the high-level template of the basic population events: a generation, as shows in Figure 1. 

As is evident, each generation comprises the updating of a population of individuals, which it is 
hopeful, will give rise to better solutions for the problem being addressed. There are three primary 
components: selection, reproduction, and replacement. 


Process Do-Generation qt pop: individual [ ]) 
variables 
breeders, newpop : Individual{ }; 
begin 
breeders <— Select-From-Population(pop) ; 
newpop «— Generate-New-Population( breeders); 


pop «— Update-Population (pop, newpop) 


Figure 1. The generational template 


As is evident, the selection uses the information provided by the fitness function, which takes 
responsibility for the population's competition. The information obtained by the temporary management 
function (fitness function) is used to evaluate the benefits of the individual in the population; then a group of 
population is selected for reproduction (reproduction is controlled by crossover and mutation) according to 
such goodwill. This selection can be made in a variety of ways. The most popular techniques are rank-based 
methods, fitness scale methods, and game-based methods (selecting individuals based on direct competition 
in a small subgroup of individuals). Replacement is closely related to this competitive aspect. 
This component ensures that the size of the group remains the same. To this end, individual in the new 
population are replaced by certain criteria for new solution. In most cases, this can be achieved by getting the 
best (guiding function) from new population (the so-called "plus" refresh strategy), or simply by capturing 
the best individual of the new pop and recording them. In pop to complete from the smallest 
("comma" strategy) [21-26]. 


3. PROPOSED MEMETIC ALGORITHM FOR MOMPS 

The Meta-algorithms are study strategies using a population-based approach where a group of 
collaborating agents and competitors participate in periods of individual enhancement to resolve 
intermittently. MA relates to a metaheuristics family whose primary theme is hybridization and is 
fundamentally interested in exploiting all accessible information about the issue under study [20]. Each local 
search algorithm investigates distinct processes for reaching individuals [22]. The suggested algorithm 
contains nine comprehensive steps: 
Step 1: Define MPS Parameters 

For MPS, the software implemented takes into consideration the following parameters (K, P, TH, 
OH, GRxp, SSK). 
Step 2: Define the Objective Function 

By using min-max normalization approach, MPS objective function can be defined as bellow: 





, EI RN. BSS oc 
Min Z = cy *—— +c * +3 * +c, * (14) 
Elmax RNMmax BSSmax OCmax 
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where, the coefficient c_1,c_2,c_3 and c_4 are used to indicate the significance of each MPS performance 
measure. El max: RNMmax, BSSmax and OCmax max are the maximum values of corresponding goals which 
are estimated from pre-processing step (warm-up period) in proposed algorithm runs. Since those values 
about objectives goals have totally different ranges, memetic algorithm is performed to place them on 
interval (0, 1). 
Step 3: Create Initial Population 

In this step, generating a population (Pop) of Np individuals. In our implementation, this population 
is randomly generated according to the uniform distributions. The pseudo code of multi (resources, products 
and periods) may be written as follows: 


for k=1:K 
for r=1:R 
if UR(k,r)4 0 
for p=1:FP 
IP=randi([0,round(GR(k,p)/BS(k,p))],nPop, 1)*BS(k,p); 
Pop=[Pop IP]; 
IP=[]; 
end 
else 
for p=1:P 
IP=zeros(nPop, 1); 
Pop=[Pop IP]; 
IP=[]; 
end 
end 
end 
end 


Step 4: Stopping Conditions 

In this algorithm, stop the procedures and return the results after a given number of iterations 
(Niter). 
Step 5: Selection of Pc and PL 

Selection is a process in which solutions are chosen according to their objective function value. 
In this algorithm, a highest rank selection is used to select population for applying crossover (Pc) and 
population for local search (PL) from Pop. In this method, we select the first individual that has greatest 
fitness (minimum Z), after sorting the population in ascending order, while the other is selected randomly. 
This method allows to improve the individuals to carry some good properties from the individual, which has 
the highest fitness. 
Step 6: Apply Arithmetic Crossover to Pc 
Output: 

The crossover produces new offspring individuals from parent individuals. 
a) Create two new individual’s crossovers by exchanging some genes of two parent individuals. 
b) Produce two new offspring by combine two linear crossovers with two parent individual vectors based 

on the following equation: 


offspringone= b* Parentone+ (1-b) * Parenttwo (1) 
offspringtwo= (1-b)*Parentone+b*Parenttwo (2) 


Where b is a random weighting factor (chosen before each crossover operation). 
Step 7: Apply Local Search Operators to PL 

The operators are the Local Search's building blocks, since they are used to travel from solution to 
solution. The choice of operators controls the richness of a solution's neighborhood. Using a large variety of 
operators increases the size of the solution space that can be reached during each iteration. This algorithm is 
using swap operator and inverse movement operator to improve the solutions of PL. 
Step 8: Evaluate Z (x) for PC and PL 

After we apply arithmetic crossover to PC and Local Search Operators to PL, we will evaluate each 
individual x in PC and PL using MPS objective function which defines in (14). 
Step 9: Select the new Population from (Pop, PC, PL) 
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The selection is carried out to obtain the next generation after combining the individuals in pop with 
PC and PL. Then elitism is ensured, and sorting is carried out on the basis of the objective function of the 
recent population. The new generation is consequently filled from the best solutions until the size of the 
population exceeds the existing size. Next, go to Step 4 


4. RESULTS AND ANALYSIS 

The production scenarios (S1, S2, S3, S4, S5, S6) are selected from Sultan [23-24], Ribas [25] and 
Supriyanto [13] to study applicability of memetic algorithm for solving MPS problem. Each scenario consists 
of (K, R and P) where K, R, P represent the numbers of production, production lines and periods 
respectively. All scenarios are tested using the software written in MATLAB 8.1 programming languages. 
The results obtained for solving MPS problem using proposed algorithm are compared with GA and SA 
algorithms as shown in Table 1. 

Master production schedules created by GA, SA and MA algorithms. Other measures, like computer 
execution time, although quite important for some types of problems, have not been considered in the present 
study. Results obtained for these scenarios, for each technique used, are shown in Table 2. The solutions of 
SA are worse than GA in all production scenarios which is the same conclusion of Viera et al. [8]. 
The Tukey’s method was used to test the results of these algorithms. In this method, the results must be 
normalized as follows: 


Normalized Ri = (Ri — min (R1: R3))/ min (R1: R3) (15) 
Where Ri is final solution of ith algorithm. These normalized finding should be evaluated using the 


methods of Tukey. The error portion of these interval plots has been presumed to be 0.05. From Figure 2, 
we notice that the normalized means of results of three algorithms have clear distinctions. 


Table 1. Production Scenarios (S1-S6) 








Production Scenario (K, R,P) Nvar(K*R*P) Source 
S1 (4,6,4) 96 Sultan et. Al. [23] 
S2 (2,9,6) 108 Sultan [24] 
S3 (4,4,7) 112 Ribas [25] 
S4 (4,4,10) 160 Supriyanto [13] 
S5 (4,4,20) 320 Supriyanto [13] 
S6 (20,4,13) 1040 Ribas [25] 





Table 2. The Computational Results of Comparison Algorithms 








PS ALG EI RNM BSS OC Z 
sl SA 1843 908 320 6.33 7.125 
sl GA 1745 854 548 4.8 6.062 
sl MA 1093 205.5 196 1.8 2.241 
s2 SA 35424 234 2517.7. 9.31 10.073 
s2 GA 35588 314 2485.7 14.33 15.101 
s2 MA 34073 123 1731.96 3.8 4.446 
S3 SA 5464.3 0 13.81 23.51 23.525 
S3 GA 5226 986 585 4.33 4.532 
S3 MA 44285 942.8 528.5 0.33 0.514 
S4 SA 1354 718.4 469 8.66 9.284 
S4 GA 1221 685 954 6.43 7.569 
S4 ME 999.1 421.5 169.4 0 0.254 
S5 SA 1165 354 75.6 4.9 5.076 
S5 GA 1201 160 172 4.13 4.383 
S5 MA 1074 0 69.5 0 0.089 
S6 SA 5562.6 1048.9 66.41 101.83 101.865 
S6 GA 5959.1 3981.5 2913 0 0.534 
S6 MA 5549 121.5 2736 0 0.430 
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The results of the comparison algorithms are thus proven according to the outputs of Tukey's 


method. Considering the execution of these three algorithms, we can clearly conclude the effectiveness of the 
proposed memetic algorithm compared to the other algorithms. 


5. 





Interval Plot of SA, GA, MA 
95% CIfor the Mean 





Data 




















Figure 2. The Confidence Interval of Three Algorithms for Final Solutions 


CONCLUSION 
This paper proposes Memetic Algorithm combines evolutionary operators with local search 


operators for solving production planning, in particular, multi-objective MPS problems and its performance 
is evaluated by different production scenarios. The MA can solve MPS ' multi-objective model effectively 
owing to its balance between exploring and exploiting a search scope. It has the ability to intelligently 
determine how much, when, and where additional capacities are required so that the inventory can be reduced 
without affecting the level of customer service. The comparison results show the efficiency and capabilities 
of the proposed algorithm in finding the optimum especially with a problem of large number of variables. 
The performance achieved is quite satisfactory and promising for solving MPS models. 
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